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Abstract 

We investigate the propagation of ultra-high energy cosmic ray nuclei (A = 1 — 56) from cosmo- 
logically distant sources through the cosmic radiation backgrounds. Various models for the injected 
composition and spectrum and of the cosmic infrared background are studied using updated photodis- 
integration cross-sections. The observational data on the spectrum and the composition of ultra-high 
energy cosmic rays are jointly consistent with a model where all of the injected primary cosmic rays 
are iron nuclei (or a mixture of heavy and light nuclei) . 
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I. INTRODUCTION 



The origin of the highest energy cosmic rays is among the most pressing questions in as- 
troparticle physics 0, 0|. The ultra-high energy cosmic ray (UHECR) spectrum has been 
measured, using both air shower and atmospheric fluorescence techniques, to energies beyond 
10 20 eV BflSi. Protons of such high energies are not expected to be deflected significantly 
by galactic or extragalactic magnetic fields jjj, but should scatter inelastically on the cosmic 
microwave background (CMB) with an attenuation length of ^ 30 Mpc, i.e. about the size of 
the local supercluster of galaxies. The observed UHECRs do not point back to any plausible 
sources within this range and have an isotropic sky distribution, so their sources, if astrophysi- 
cal, must be cosmologically distant. If so the spectrum should be suppressed at energies above 
Eqzk ~ 6x 10 19 eV (the "Greisen-Zatsepin-Kuzmin cutoff"), if the primaries are indeed protons 
[H, [[|. Presently, there are conflicting claims concerning the existence of this spectral feature. 
Data from AGASA (Akeno Giant Air Shower Array) show no indication of any suppression 0, 
while data from the HiRes (High Resolution Fly's Eye) air fluorescence detectors are consistent 
with the expected cutoff p. The Pierre Auger Observatory, which employs both techniques, 
has accumulated an exposure comparable to AGASA and about half that of HiRes. Its results 
are presently consistent with either possibility and forthcoming data should be able to 
settle the issue definitively. The AGASA observation that the spectrum continues smoothly 
beyond the GZK cutoff has motivated many suggestions for the origin of the trans-GZK events 
such as decaying superheavy dark matter particles in the Galactic halo 2 and ultra-high 
energy neutrinos creating 'Z-bursts' in local interactions with the cosmic neutrino background 
[l3L Qj| . It has also been suggested that the primaries may be neutrinos with large interaction 
cross-sections at ultra-hi gh e nergies [lil or perhaps new stable strongly interacting particles 
heavier than the proton |l7| . 

Rather less exotic ways of evading the GZK cutoff are also possible. In particular, if UHECRs 
are not protons but consist of heavy nuclei [3] > then the spectrum will be altered from the usual 
expectation. Although evidence for the presence of heavy nuclei in UHECRs has been around 
for quite some time 2, III 21], the information available on UHECR composition is still rather 
imprecis e 1221 . t he most reliable result being that the primaries are not mostly photons but 
hadrons |23l. 124 I25L liil] . Despite the widespread impression to the contrary, however, there is 
little reason to believe that these particles are protons rather than heavier nuclei. 

Determination of the mass composition of UHECRs has been attempted using a variety of 
variables such as the rate of change with primary energy of the depth of shower maximum 
A max , the degree of fluctuation in the depth of the shower maximum, the number of muons 
which reach the Earth's surface, the lateral density profile with respect to the shower axis, the 
width of cosmic ray shower disks, etc. Such measurements must be compared to Monte Carlo 
simulations, typically conducted with programs such as QGSJET, SIBYLL and DPMJET |27| . 
The mass composition inferred is quite sensitive to the simulation program used. Furthermore, 
data from different experiments can yield rather different results regarding cosmic ray composi- 
tion even when the same shower simulation is used. For example, at energies near 10 18 eV, the 
mass composition determined from the lateral density profile of Volcano Ranch and Haverah 
Park data are in disagreement at > 2a significance [28(. The situation does not improve much 
at higher energies. Above 10 19 eV, only the Fly's Eye/HiRes and AGASA experiments have 
reported any information pertaining to mass composition, and these results extend to only 
about 2.5 x 10 19 eV — in this energy range the data suggest a predominantly light composition 
although the uncertainties are large jjj ]. 
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So far no individual sources of UHECRs have been identified. This can be understood if 
there are a very large number of faint sources or, alternatively, if there are large-scale magnetic 
fields which deflect UHECRs sufficiently so as to to conceal their origin. As mentioned above, 
protons at trans-GZK energies are not expected to be significantly deflected (but see ref. 1 301) . 
Heavy nuclei, by contrast, have greater electric charge and may thus have their arrival directions 
isotropised, even at such high energies. 

There is also a theoretical prejudice to expect the UHECR spectrum to be dominated by 
heavy nuclei, since the maximum energy to which particles can be accelerated scales with the 
electric charge. Thus astrophysical candidates for UHECR accelerators such as active galactic 
nuclei and gamma-ray bursts, which barely meet the "Hillas criterion" j3j| for containment and 
acceleration of 10 20 eV protons, can in principle easily accelerate heavy nuclei to such energies. 

Given these arguments, and the inconclusive nature of the experimental data regarding 
UHECR composition, it is rather surprising that most studies in this field have usually as- 
sumed that UHECRs are protons. By contrast , relatively little attention has been paid to the 
propagation of ultra-hi gh e ner g y nucle i |32L I33L 34l | although there has been a resurgence of 



interest in recent years |35|, |36|, |37|, |38|, |39|, |4£], |4l|. We revisit this problem, using up-to-date 



nuclear physics data, and considering a wide range of nuclei as possible UHECR primaries. We 
also explore different models for the cosmic infrared background (CIB), as well as the effects of 
possible intergalactic magnetic fields. We pay particular attention to the relationship between 
the composition of UHECRs at Earth and that injected at source. 

We begin by considering the intergalactic propagation of ultra-high energy (UHE) protons 
(Section [nj. In Sections llll| I1VI and IV! we study the photodisintegration of cosmic ray nu- 
clei, including the effect of uncertainties in the cross-sections and the CIB model used. We 
then proceed to discuss the composition of UHECRs observed at Earth (Section |VTJ), and the 
role of magnetic fields in the propagation of cosmic ray nuclei (Section IVII|) . We propose 
some representative models of UHECR primary composition and discuss X max measurements 
in Sections IVIIII and IIXI Our conclusions are presented in Section 1X1 



II. PROPAGATION OF ULTRA-HIGH ENERGY PROTONS 

Before discussing the propagation of intermediate mass and heavy nuclei, we first recapitulate 
the physics governing the propagation of UHE protons. 

Over cosmological distances, the dominant processes effecting UHE proton propagation are 
interactions with the CMB producing pions or electron-positron pairs. Pair production (p + 
7cmb p + e + + e~) occurs sufficiently frequently that it can be treated as a continuous energy 
loss process |42j ; as shown in Figure Q this effect dominates over the energy loss due to the 
Hubble expansion for proton energies ~ 10 18 — 10 21 eV and peaks at ~ 3 x 10 19 eV. 

Pion production is not as simple to implement. Individual occurrences of the processes 
P+7cmb — > p+7T° and p+7cmb — ¥ n+ir + cause the primary proton to lose a considerable fraction 
of its energy. Therefore these cannot be treated as a continuous processes and it becomes 
necessary to use Monte Carlo techniques 43||. Furthermore, if enough energy is available, 



multi-pion production can be important; these non-perturbative processes take place through 
the near- resonance exchange of the 1.232 GeV A + -hadron 0|. The associated energy loss 
lengths are shown in Figure ^ 

It is also possible for UHE protons to produce pions through interactions with the cosmic 
infrared background (CIB) as shown in the right frame of Figure ^ We show results for three 
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FIG. 1: Energy loss lengths for UHE protons propagating through the universe. The left frame 
shows the results for the processes p + 7cmb — ► P + e + + e~ ("e + e~ creation"), p + 7cmb — ► P + 
P + 7cmb — ► n + 7r + and multi-pion production ('V production"), as well as due to the Hubble 
expansion ( "redshift" ) . In the right frame, only pion production processes are considered and the 
effects of the cosmic infrared background (CIB) are shown in addition to the CMB (the three CIB 
models are discussed in Section ITV)l . 



models of the CIB which are discussed further in Section IIVI Although the rates for these 
interactions are sub-dominant in comparison to energy losses from pair production, they can 
be important in determining the spectrum of "cosmogenic" neutrinos produced in UHE proton 
propagation [45| . 

For the injection spectrum, we follow previous work in adopting a simple power law with a 
cutoff: 

^ ex E p < (£ max /26). (1) 

Later in this paper, the quantity E max will be used to describe the maximum energy to which 
iron nuclei can be accelerated by astrophysical sources. Since this maximum energy scales with 
the electric charge of the accelerated particle, protons have a cutoff energy 26 times smaller. 
We also assume a constant comoving density of sources: 

The possible effects of source number evolution has also been considered, typically by adopting 
a distribution proportional to (1 + z) 4 . We find however that this has a negligible impact as 
UHECRs typically propagate over distances of only 10-100 Mpc (i.e. z 1). In Figure we 
show the UHECR spectrum expected at Earth for proton primaries — the "GZK cutoff" is 
clearly seen. It is also evident that none of these models fit the data well. 



III. PHOTODISINTEGRATION OF ULTRA-HIGH ENERGY NUCLEI 

The propagation of UHECRs is quite different in the case of heavy nuclei. These undergo 
photodisintegration in scattering off the CMB and/or CIB at a rate: 

A 2 mlc 2 r°°den(e) r 2Ee / Am P c 



RA,z,i p ,i n - 2 ^ J — ~2 — L de'e'(TA,z,i p ,i n (t'), (3) 
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FIG. 2: The UHECR spectrum at Earth for purely proton primaries injected by homogeneously 
distributed extragalactic sources with power law spectral indices of a = 2.0, 2.4 and 2.8 up to a 
maximum energy of (S max /26]_ = (10 22 eV/26) ~ 10 20 ' 5 eV. The overall flux has in each case been 
normalized to the Auger data ilOj (the x 2 /d.o.f. for the fits are 9.97, 4.36, and 1.49 respectively). The 
effects of magnetic fields are not included. 



where A and Z are the atomic number and charge of the nucleus, i p and i n are the numbers 
of protons and neutrons broken off from a nucleus in the interaction, n(e) is the density of 
background photons of energy e, and crA,z,i p ,i n ( e> ) is the appropriate cross-section. (Note that 
e is measured here in the laboratory rest frame, but in subsequent equations it refers to the 
energy measured in the rest frame of the nucleus.) 

The cross-sections for photodisintegration have often been modelled using the parameteri- 



zation of Stecker and collaborators 134, 46, 4 
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where £i, C> e p,i an d Aj are parameters whose values are obtained by fitting to nuclear data. 
Here, i is the total number of nucleons broken off from the nucleus in the interaction and the 
integrated cross-section is 
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while the function Wi is given by 
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Here Q+(x) and 9_(a;) are the Heaviside step functions, e x = 30 MeV, e max = 150 MeV, and 
the threshold energy for a given process is in most cases e t h r ~ i x 10 MeV (values are tabulated 
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FIG. 3: Energy loss lengths for photodisintegration of oxygen (left) and iron (right) nuclei on the 
CMB and CIB for both Gaussian and Lorentzian parameterizations of the cross-section. The CIB has 
been modelled according to Malkan & Stecker |48| (see Section II V|) . 



in Ref. 34]). These cross-sections are dominated by the giant dipole resonance which peaks 
in the energy range ~ 10 — 30 MeV; at higher energies, quasi-deuteron emission is the main 
process. 

This model incorporates two major approximations. First, a simple Gaussian form is as- 
sumed for the photodisintegration cross-section, cut off abruptly below the theoretical reaction 
threshold. Second, a fixed choice of Z is made for each given atomic number, A, thus neglecting 
the many possible isotopes which may be generated in such interactions. 

Although this parametrization does a fairly good job of reproducing the measured cross- 
sections, possible improvements have been proposed. In particular, a generalized Lorentzian 
model, which can be fitted to the available nuclear data, has been used to parametrize a wide 
range of photodisintegration cross- sect ions j3?J: 
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resonance, given by 
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is the position of the giant dipole resonance and TA,z,i p , 
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In some cases, the Lorentzian form can be quite different from the Gaussian parameteri- 
zation. In particular, the latter often overestimates the width of the giant dipole resonance. 
Despite these differences we find, as shown in Figure EJ that the energy loss rates due to pho- 
todisintegration change little whether we use the Gaussian model or the Lorentzian model 
j3?| for nuclei between A = 11 and 56 — henceforth we adopt the Lorentzian model. (Below 
this mass range, we use the Gaussian parameterization.) 



IV. THE COSMIC INFRARED BACKGROUND 



Since photodisintegration processes occur most efficiently when the Lorentz-boosted target 
photon can excite the giant dipole resonance at ~ 10 MeV, an UHE nucleus with an energy 
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FIG. 4: Representative models of the CIB spectrum from Malkan & Stecker |48l ]. Aharonian a£. 
and Franceschini et <rf. 0. Data shown are from DIRBE 51], ISO 0, HST 0, El and 
IRTS 0. 



of 10 19 eV (i.e. a Lorentz factor of ~ 10 s — 10 9 ) needs to collide with a ~ 0.01 — 0.1 eV 
background photon in order to most efficiently undergo photodisintegration. Photons of this 
energy, corresponding to wavelengths of A ~ 10 — 100 fim, are present in the CIB rather than 
in the CMB. 

The CIB is an expected relic of the cosmological structure formation processes The 
assembly of baryonic matter into stars and galaxies and the subsequent evolution of such 
systems is accompanied by the release of radiant energy; cosmic expansion and the absorption 
of short wavelength radiation by dust and re-emission at longer wavelengths shifts a significant 
part of this radiant energy into the infrared: A ~ 1-1000 /im. Thus the CIB spectrum depends 
on the luminosity and evolution of the sources and distribution of dust from which it is scattered. 

Direct measurement of the CIB has been performed by several satellites, in particular the 
COsmic Background Explorer (COBE) and the InfraRed Telescope inSpace (IRTS). DIRBE, an 
instrument aboard the COBE satellite (52|, has provided measurements in the wavelength range 
1.25-240 /im. FIRAS, another instrument aboard COBE, covered the range 25 /im - 1 mm. 
The ISO satellite carried two instruments which were employed in the indirect measurement of 
the CIB: ISOCAM at 7 and 15 /im and ISOPHOT at 170 fim [HJ. 

In addition to these measurements, telescopes such as the Hubble Space Telescope's (HST) 
wide field planetary camera, combined with spectrophotometry from the duPont telescope and 
the HST Faint Object Spectrograph, were used to measure the CIB at 0.3, 0.55 and 0.8 /xm. 
Galaxy counts made using the HST Northern and Southern Deep Fields, between 0.36 and 2.2 
/im, supplemented with shallower ground based observations, have also been used to set lower 
limits on the CIB 0, H • 

We have considered three forms for the CIB spectrum which bracket the range of pos sibilities, 
The first is based on a compilation of the direct observations by Aharonian et al |49( which is 
the uppermost curve at 1 fim in FigureEJ The second, from Franceschini et al. |5CI |. corresponds 
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FIG. 5: Energy loss lengths due to photodisintegration for helium, oxygen and iron nuclei for three 
models of the CIB spectrum. The Lorentzian model (3^ for photodisintegration cross-sections, been 
used. 



roughly to the lower limit set by galaxy counts. The third is the empirical model by Malkan 
& Stecker |48[ which lies between the other two curves (although at 100 fim it goes below the 
Franceschini et al. model). The observational data are also shown for comparison. 

In Figure El we show the effect of our choice of the CIB spectrum on the energy loss rates of 
cosmic ray nuclei through photodisintegration. It is clear that the choice makes a substantial 
difference only for particularly heavy nuclei. Henceforth, unless otherwise stated, we will adopt 
the Malkan & Stecker model [48[ for the CIB spectrum. We note that this spectrum is consistent 
with recent observations by HESS of TeV 7-rays from distant blazars which imply restrictive 
upper bounds on the CIB and suggest that the direct observations in the ~ 1 — 10 /im range 
may well have been contaminated by zodiacal light EF 



V. PROPAGATION OF ULTRA-HIGH ENERGY NUCLEI 

Now we can study the intergalactic propagation of UHE nuclei and determine the UHECR 
spectrum at Earth for various types of nuclei injected at source. In Figure |H1 we plot the 
energy loss length due to photodisintegration for several nuclei species. It is seen that there are 
significant variations, for example, carbon nuclei are relatively robust to photodisintegration 
with an energy loss length of ~ 50 Mpc at 10 20 eV and several thousand Mpc at 3 x 10 19 
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eV, while beryllium nuclei are highly fragile, with an energy loss length smaller by a factor of 
~ 10 — 100. Very heavy nuclei are generally quite stable up to energies ~ 10 20 eV. 

In practice, a heavy nucleus would undergo many photodisintegration reactions, cascading 
down in atomic number and charge. Thus over a given trajectory, an injected particle will 
take on the identity of many species of nuclei, generating UHE protons, neutrons and alpha 
particles along the way, each of which will continue to propagate and contribute to the UHECR 
spectrum. In Figured we show the spectrum observed at Earth for various species of injected 
nuclei, assuming an injection spectrum: 

^ oc E~ a } E < (£ max x Z 1 26), (9) 
where Z is the charge of the nuclei species under consideration. 



VI. THE COMPOSITION AT EARTH 

If intermediate mass or heavy nuclei are injected in a distant cosmic ray accelerator, these 
particles will gradually disintegrate into lighter nuclei and nucleons as they propagate through 
intergalactic space. Depending on the distance to the sources, the cosmic ray composition 
observed at Earth may be quite different from that at injection. 
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FIG. 7: The spectrum of UHECRs observed at Earth for a range of injected heavy nuclei with power- 
law spectral index a = 2.4 or 2.0 and E mm = 10 22 eV. The overall flux has in each case been 
normalized to the Auger data ^3). The Malkan & Stecker CIB model (3] and the Lorentzian model 
[37] | for photodisintegration cross-sections have been used. The effects of magnetic fields have not been 
included. 10 



The observed composition at Earth has a distinctive dependence on the energy as can be 
seen in Figure |HJ In the interesting energy range ~ 3 x 10 19 — 10 20 eV, i.e. where the GZK 
suppression is expected for proton primaries, photodisintegration is most effective. For sources 
injecting intermediate mass nuclei, the effective mass number at Earth reaches a well-defined 
minimum, probably indistinguishable from proton primaries. This minimum is less pronounced 
for injection of very heavy nuclei and the effective composition at Earth is distinctly heavier 
than protons. 

These issues are of critical importance, observationally speaking. The injected composition 
of the UHECR spectrum is not directly accessible experimentally, and can only be reconstructed 
from the composition observed at Earth. As stated earlier, the present observational status is 
rather uncertain. Future data will hopefully reach a level of quality which makes it possible to 
reliably infer the approximate composition at Earth. With such data, general trends such as 
those seen in Figure |H] would aid in estimating the composition of cosmic rays at the sources. 



VII. EFFECTS OF INTERGALACTIC MAGNETIC FIELDS 

So far in this study, we have neglected the effects of magnetic fields on the propagation of 
UHECRs. For protons or nuclear primaries, however, such effects can play an important role 
in determining the cosmic ray spectrum. The importance of these effects depend, of course, on 
the strength of the extragalactic magnetic fields which is currently a subject of some debate 
with contrary conclusions drawn in Ref. 0| and in Refs. |30l l38l|. 

A charged particle moving through a uniform magnetic field undergoes an angular deflection 
upon traversing a distance, L co h, of a = L^/Rl, where Rl is the Larmor radius of the particle. 
Therefore a particle traversing a distance, L, through a series of L/L co h randomly orientated 
uniform magnetic field regions of length L co h, suffers an overall angular deflection given by 

where L co ^ is the representative coherence length of the extragalactic magnetic fields, B is their 
representative magnitude and Z is the electric charge of the cosmic rays. Such deflections result 
in an increase in the effective distance to a cosmic ray source given by: 

L y 1 ! 2 V E ) UOMpcAlMpcAlnG/ V267 v ; 

Thus for protons or light nuclei, nano-Gauss magnetic fields have little impact for the high 
energies considered here. This is not true for heavy nuclei, e.g. for iron nuclei propagating 
through nG-scale magnetic fields, the effective distance to a source 50 Mpc away is increased 
by ~ 30% at 10 20 eV (alternatively, the energy loss length is reduced by about ~ 30%). Since 
this effect scales with the inverse square of the cosmic ray energy, such (plausible strength) 
magnetic fields would have a dramatic effect on the propagation of lower energy heavy nuclei. 

In FigureElwe show the effects of such extragalactic magnetic fields on the UHECR spectrum. 
For oxygen primaries, the effects are small, only becoming of any consequence at energies below 
a few times 10 19 eV. However the effects are more prominent for iron primaries. 

Some words of caution are called for at this point. The effects of nG-scale magnetic fields 
appear to set in at an energy of roughly 5 x 10 19 eV for oxygen, whose primaries can arrive 
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FIG. 8: The mean atomic mass of cosmic rays arriving at Earth for a range of injected heavy nuclei 
with power-law spectral index a = 2.4 or 2.0 and £ max = 10 22 eV. The Malkan & Stecker CIB model 
ji^l and the Lorentzian model for photodisintegration cross-sections have been used. The effects 
of magnetic fields have not been included. 
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FIG. 9: The effects of nano-Gauss extragalactic magnetic fields on the UHECR spectrum for oxygen 
and iron primaries with power-law spectral index a = 2.4 and £' max =10 22 eV, assuming -L C oh ~ 1 Mpc. 
The overall flux has in each case been normalized to the Auger data 10]. The Malkan & Stecker CIB 
model [3 and the Lorentzian model [s^J for photodisintegration cross-sections have been used. 



from approximately 300 Mpc at this energy, corresponding to an effective length of roughly 
L e s/L ~ 1 + 0.7 x (B/nG) 2 . This means that the small angle treatment is valid for magnetic 
fields of strength upto ~ 0.3 nG. On the other hand, an iron nucleus with an energy of 5 x 10 19 
eV could have traveled hundreds or even thousands of Mpc before its arrival at Earth, so would 
have been deflected by an angle of 9 ~ 130° x (L/100 Mpc) 5 x (B/nG). Hence the result in 
the right frame of Figure |H1 is not reliable at low energies; to properly take into account such 
strong deflections, a numerical simulation of diffusion including an appropriate description of 
the magnetic field structure in the local supercluster is required. Such a treatment is however 
beyond the scope of the present study. 



VIII. MIXED ULTRA-HIGH ENERGY COSMIC RAY COMPOSITION 

We now study the possibility that a mixture of protons, helium, oxygen, and iron nuclei 
are injected at source in roughly the same proportions as those observed in low-energy galactic 
cosmic rays In particularl we will consider a mixture in the ratio of H : He : O : Fe 

= 1 : 0.85 : 0.06 : 0.02 for a power-law spectral index a = 2.0, and H : He : O : Fe = 
1 : 1.48 : 0.19 : 0.09 for the case of a = 2.4. 

The injection spectrum is assumed to have a cutoff modelled as: 

E \/26- 



— ocE- ex P 



(12) 



max - 

with the cutoff energy in the range: 10 20,5 eV < E max < 10 22 eV. Once again, the choice of 
the spectrum is motivated by the possibility that cosmic ray sources can accelerate charged 
particles to a maximum energy proportional to the charge. This leads to the expectation that 
such sources will preferentially accelerate heavy nuclei to the highest energies, if such particles 
are indeed able to survive the radiation fields present in the sources. 

In the left frames of Figure ITUl we plot the propagated spectrum calculated using this mixed 
composition at source, compared to that for pure proton injection. It is seen that the two cases 
are in fact hard to distinguish observationally. 
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FIG. 10: The spectrum (left) and depth of shower maximum, X max , (right) for a mixture of protons, 
helium, oxygen, and iron injected at source with proportions following that of the Galactic cosmic 
rays. The solid, dashed, and dotted lines are the predicted values of X max for, respectively, a pure iron, 
pure proton, and the mixed composition, obtained using three simulation programmes — DPMJET, 
SIBYLL and QGSJET. The overall flux has in each case been normalized to the Auger data [lo| . 
The X max data are from the Fly's Eye 0|, HiRes Epl a nd Yakutsk [2(J experiments. The Malkan 
Sz Stecker CIB model and the Lorentzian model|3^| for photodisintegration cross-sections have 
been used. 



In the right frames of Figure EH we plot the average value of X max as a function of the 
energy for the same mixed composition at source. To illustrate the uncertainties in relating 
Xmax to the UHECR composition, we 



lave shown the results of three air shower simulation 
62j, and QGSJET [61]. Also shown in Figure EH are 
IS |, HiRes |2^ |. and Yakutsk J2(j experiments. (The 



programmes: DPMJET |61|, SIBYLL 
the measurements from the Fly's Eye 
Xmax values quoted in the Auger analysis which sets an upper bound to the photon content in 



UHECRs |26j are explicitly stated to be inappropriate for elongation rate studies.) The data 
favour a composition heavier than protons, regardless of which simulation program is adopted. 



IX. PREDICTIONS FOR X max MEASUREMENTS 

Information on the UHECR composition at Earth can be obtained from studies of cosmic 
ray shower development, in particular the lateral distribution, muon content, and X max value 
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22l | . In this Section, we will focus on the parameter X max , which is the atmospheric depth at 



which the number of particles in the shower is largest. Measurements of X max are made by 
fluorescence detector experiments which measure the UV light emitted by the excited atoms 
(mainly nitrogen) in the air shower. 

In the left frames of Figure ITU we show the UHECR spectrum and values of X max calculated 
for iron, oxygen and helium nuclei injected at source, with a power-law spectral slope of a = 2.4 
and -E max — 10 20 ' 5 eV. With this rather low choice of the cutoff energy, heavy nuclei dominate 
the spectrum above ~ 10 19 eV — iron and oxygen nuclei break up rather modestly during 
propagation, so the composition at Earth is similar to that at source. In the right frames of 
Figure ^TJ we repeat the calculations for the case of a = 2.0 and find similar results.. 

For the case of a low cutoff energy, we expect a transition from a heavy composition at high 
energies to a much lighter composition at low energies, largely independently of the ratio of 
species present in the source environment. From the spectra shown in Figure ^JJ we would 
expect a nearly all- iron composition at the highest energies, and a nearly all proton/helium 
composition below ~ 10 19 eV, with the value of X max changing over this range accordingly. 
The current data on X max is too scattered to establish whether such a transition is present. 

This conclusion can be modified if the accelerators of UHECRs operate up to a higher cutoff 
energy so that the heavy nuclei injected at source are substantially photodisintegrated by the 
time they reach Earth. In Figure we show that for the case £ max = 10 22 eV, it is not at 
all clear which species will dominant the observed UHECRs — the data is equally well fitted 
by an all-proton, all-iron or mixed composition at source. What is interesting is that the X max 
data is reasonable well fitted when only oxygen nuclei are accelerated by the sources. Even 
pure iron injection at source matches the data for an injection spectrum with a = 2. 



X. DISCUSSION 

We have studied the intergalactic propagation of a variety of heavy and intermediate mass 
cosmic ray nuclei at ultrahigh energies. Adopting different models for the cosmic infrared 
background and for the photodisintegration cross-sections has little effect on the propagated 
energy spectrum and composition at Earth. Of more significance is the choice of the source 
spectrum and the effect of intergalactic magnetic fields. Our main aim was to determine 
the relationship between cosmic ray composition at source and at Earth after the effects of 
propagation are taken into account. Somewhat surprisingly, extant data on the composition of 
UHECRs is consistent with the injection of even pure iron nuclei by the sources. 

As the Pierre Auger Observatory continues to accumulate more exposure to UHECRs, a 
definitive resolution is expected soon of whether there is indeed a GZK cutoff in the spec- 
trum. By combining the energy spectrum obtained using the surface detectors with X max 
measurements by the fluorescence telescopes, information can then be extracted on the sources 
of UHECRs using the results presented in the present paper. Observations of other types of 
messengers associated with the highest energy cosmic rays, such as the cosmogenic neutrino 
flux, will also help to determine whether the ultra-high energy cosmic rays are protons or heavy 
nuclei jH^] amd bring us closer to answering the long standing mystery of their origin. 
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FIG. 11: The expected spectrum and X max for helium, oxygen, and iron nuclei injected at source for 
E max =10 20 ' 5 eV with a = 2.4 and 2.0. The overall flux has in each case been normalized to the Auger 
data 10]. The X max data are from the Fly's Eye 0], HiR.e s and Yakutsk ji3] experiments. The 
Malkan & Stecker CIB model [3] and the Lorentzian model |37i] for photodisintegration cross-sections 
have been used. 16 
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